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Abstract 



^| The charm-loop contribution to B — > X s ^ is found to be numerically dominant and very 
i-^h ■ stable under logarithmic QCD corrections. The strong enhancement of the branching ratio 

^ ■ by QCD logarithms is mainly due to the 6-quark mass evolution in the top-quark sector. 
^ . These observations allow us to achieve better control over residual scale-dependence at the 

^ \ next-to-leading order. Furthermore, we observe that the sensitivity of the matrix element 
(X s ^\(sc)v-A{cb)v-A\b) to m c /mb is the source of a sizeable uncertainty that has not been 
properly taken into account in previous analyses. Replacing m^ ole /ml° le in this matrix element 
by the more appropriate m^ s (/z)/m£ ole with /i G [m c ,mb] causes an 11% enhancement of the 
SM prediction for BR[i? — > X s j]. For E 1 > 1.6 GeV in the 5-meson rest frame, we find 
BR[£? — > X s 7]p>i.6 GeV = (3.60 ±0.30) x 10~ 4 . The difference between our result and the cur- 
rent experimental world average is consistent with zero at the level of la. We also discuss the 
implementation of new physics effects in our calculation. The Two-Higgs-Doublet-Model II 
with a charged Higgs boson lighter than 350 GeV is found to be strongly disfavoured. 



* Comparison with experiment is updated in the present hep-ph version. In particular, the published results of 
CLEO U are taken into account. 



1 Introduction 



Strong constraints on new physics from B — > [[j]-|3fl crucially depend on theoretical un- 
certainties in the Standard Model prediction for this decay.Q This becomes more and more 
transparent with the progress in experimental accuracy. The current experimental results are 
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BR[B^X sl \ = 3.21±0.43 st at±0.27 ; 

BK[B^X sl \ = [3.36 ± 0.53 stat ± 0.42 sys 
BR[6^s 7 ] = (3.11±0.80 stat ±0.72 sys 

The weighted average for BR 7 =BR[i? — > X s / ~f] is therefore?] 

BR!* P = (3.23 ± 0.42) x 10~ 4 
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with an error of around 13%. Improved results from CLEO and BELLE, as well as new results 
from BABAR, are expected soon. Thus, all efforts should be made, on the theoretical side, to 
reduce the uncertainty significantly below 10%. 

There are three sources of uncertainties in the theoretical prediction: parametric, non- 
perturbative and perturbative. The most important parametric ones in our analysis are due to 
a s (Mz) and the HQET parameter Ai. The latter parameter occurs in the calculation of the 
phase-space factor of the semileptonic branching ratio that is conventionally used for normal- 
ization. 

All the non-perturbative effects, except for the calculated A 2 /m 2 cb corrections [|7|,§]> sum U P 
to the non-perturbative uncertainty. No satisfactory quantitative estimate of those effects is 
available, but they are believed to be well below 10% when the photon energy cut-off is above 
1 GeV and below 1.9 GeV jnj in the 5-meson rest frame. The cut-offs imposed most 
recently by CLEO and BELLE are 2.0 and 2.1 GeV, respectively, in the T(4S) frame. 

As far as the perturbative uncertainties are concerned, they were dramatically reduced 4 
years ago, after the completion of NLO QCD calculations ||Tl~HT^]. A further improvement 
came from electroweak corrections [|i~0|. |l7|-|i~9|l . The unknown NNLO QCD corrections were 



estimated at the level of ±7% in ref. flOfl , where it was pointed out that certain scale-dependence 
cancellations at NLO seem accidental. 

However, one source of perturbative uncertainty was not properly taken into account in the 
previous analyses |l0|,|l^ , [19| , p0l . It is related to the question of the definitions of m c and nib that 
should be used in the matrix element (P2) = (X s j\(sc)v-A(cb)v-A\b) ■ This matrix element is 
non-vanishing at two loops only, so the renormalization scheme for m c and m\, is an NNLO 



1 B denotes cither B° or B~ , while X s stands for S = — 1 hadronic states containing no charmed particles. 

2 Statistical errors in the ALEPH measurement of b — * S7 are much larger than expected differences among 
weak radiative branching ratios of the included 6-hadrons. 
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issue (unless In ^ is treated as a large logarithm). However, this problem is numerically 
very important because of the sensitivity of (P2) to m c /mb. Changing m c /rrib in (P 2 ) from 
0.29 ± 0.02 to 0.22 ± 0.04, i.e. from mv° lc /ml° le to mp(/i)/m£ ole (with // G [m c , m h \) causes 
an increase of BR 7 by around 11%. What matters for B — > Xg'j is mainly the real part of (P2), 
where the charm quarks are usually off-shell, with a momentum scale set by m^ ole (or some 
sizeable fraction of it). Therefore, as we shall argue in appendix D, the choice of m^ IS (/i)/m^ ole 
with /i G [m c ,mfe] seems more reasonable than mP° le /m^ ole . 

Once m^ IS (/i)/mjj' olc with fi G [m c ,mb\ is used in (P2), the uncertainty in BR 7 significantly 
increases. This is due in part to a strong scale-dependence of m c (/z). Moreover, in all the 
previous analyses, the m c -dependence of (P2) cancelled partially against that of the semileptonic 
decay rate. Once the different nature of the charm mass in the two cases is appreciated, the 
cancellation no longer takes place. 

In the present paper, we perform a reanalysis of B — > X s 7, taking the above problems with 
m c /mb into account. At the same time, we make several improvements in the calculation, which 
allows us to maintain the theoretical uncertainty at the level of around ~ 10%. In particular, 
a careful calculation of the semileptonic phase-space factor is performed, after expressing it in 
terms of an observable for which the NNLO expressions are known. Moreover, good control 
over the behaviour of QCD perturbation series in B — > X s 7 is achieved by splitting the charm- 
and top-quark-loop contributions to the decay amplitude. The overall factor of rrib is frozen 
at the electroweak scale in the top contribution to the effective vertex rrib^SLcr^b^F^. All 
the remaining factors of rrib are expressed in terms of the bottom "IS mass". As argued in 
ref . [21] , expressing the kinematical factors of nib in inclusive S-meson decay rates in terms of 
the IS mass improves the behaviour of QCD perturbation series with respect to what would 
be obtained using m^ s (mfc) or ml° le . When such an approach is used, no sizeable accidental 
cancellations of scale-dependence in the NLO expressions for BR 7 are observed any more. 

Splitting the charm and top contributions to the amplitude allows us to better understand 
the origin of the well-known factor of ~3 enhancement of BR 7 by QCD logarithms. When 
the splitting is performed at LO, the charm contribution is found to be extremely stable under 
QCD renormalization group evolution. The logarithmic enhancement of the branching ratio 
appears to be almost entirely due to the top-quark sector. It can be attributed to the large 
anomalous dimension of the 6-quark mass. 

Our paper is organized as follows. In section [2|, splitting of the charm and top contributions 
is performed at the leading- logarithmic level, and rrib(fi) in the top sector is shown to be the 
main source of large QCD logarithms. In section the NLO formulae for BR 7 are written in 
such a way that the NNLO uncertainties can be conveniently controlled. Section £| is devoted 
to the numerical analysis, our final result being given in eq. ( |4.14| ). In section |5|, constraints on 
new physics are briefly discussed. Section || contains our conclusions. 
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Details on specific points of our calculation as well as the formulae that are necessary to make 
our analysis self-contained are collected in the appendices. In appendix A, our input parameters 
are collected. Appendix B is devoted to a determination of the mass ratio m^' IS '(m t ) jm\ . 
In appendix C, we calculate the semileptonic phase-space factor. Appendix D contains the 
determination of the ratio m c /mi, that should be used in (P2). We also give there the analytical 
dependence of this matrix element on m c /mb [IT]. In appendix E, the relevant bremsstrahlung 
formulae are summarized. 



2 Leading-order considerations 

The resummation of large QCD logarithms in B decays usually begins by decoupling the heavy 
electroweak bosons and the top quark. In the resulting effective theory, flavour- changing in- 
teractions are present only in operators Pj of dimension > 4 . Their Wilson coefficients Ci(fi) 
evolve according to the Renormalization Group Equations (RGEs) from the matching scale Ho 
down to the scale fit, where the matrix elements of Pj are evaluated. 
In the case of B — > X s 7, the relevant operators Pj readQ 

Pi = (s Ll ,T a c L )(c L ^T a b L ), 
P2 = (s L J^c L )(c L ^b L ), 

Ps = (*L7A)E(<77^), 
q 

P4 = (s L ^T a b L )Y,(<irT a q), 
q 

P5 = is L -i^ v -i P b L )Y J i.qi ll i v i p q), 

q 



P 7 = -^—m b ^)(s L a^b R )F, u , 



Ps = ^ 2 mM{sL^ u T% R )Gl v . (2.1) 

In the leading logarithmic approximation, the b — > 37 amplitude is proportional to the (effec- 
tive) Wilson coefficient of the operator P7. The well-known expression for this coefficient 
reads 

c? )cS (^) = v%c?>M + 1 - ,1) cf („ ) + E (2.2) 

6 1=1 



3 The CKM-sup pres sed (su)(ub) analogues of Pi and P2 are present in the effective theory, too. Their effects 
are included in eq. ( |2.2| ) and in our NLO analysis in the following sections. 
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where r\ = a s (/io) / o: s (/i&) and 
hi = 



626126 
272277 



56281 
51730 



i -0.6494 -0.0380 -0.0185 -0.0057 



(2.3) 



The powers dj are given in table p], in section |3|. The coefficients (/io) and Cg (//q) are 
found from the one-loop electroweak diagrams presented in fig. |1| It is sufficient to calculate 
the 1PI diagrams only. 






G 



— V- 



w, c, t s o G s b u,c,t s 

Figure 1: One-loop 1PI diagrams for b — > S'j in the SM. There is no W ± G T ^f coupling in the 
background-field gauge. 

Contributions from different internal quark flavours in those diagrams can be separately 
matched onto gauge-invariant operators, even when the calculation is performed off-shell.^ For 
the operator Pi and its gluonic analogue Ps, each quark flavour yields a UV- finite contribution 
that depends neither on the renormalization scheme nor on the gauge-fixing parameter. 

When such a separation of flavours is made, and the CKM-suppressed w-quark contribution 



is neglected, eq. (|2.2|) can be written as 

cf )cE ^ b )=x c + x u 

where the charm-quark contribution is given by 



23 16 8/14 16 
7)23 — — ( 77 23 — r)23 

3<S ' Q V ' ' 



36 ' 9 
and the top-quark one reads 



m t 

M. 



16 
fj 23 



J- n 



i=l 



m t 



(2.4) 



(2.5) 



(2.6) 



where 



Al(x) 



3x 3 +2x 2 1 , -22x 3 +153x 2 -159a:+46 

2(x-l) 4 111 36(2-1)3 : 



2(*-l) 

74 lnx + 



2(x-l) 4 1 12(x-l) 3 

The first two terms in X c Q2.5| ) are obtained from eq. 



(2.7) 

by the following replacements: 



4 In an off-shell calculation, use of the background-field gauge is necessary to ensure the absence of gauge- 
non-invariant operators. 
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C\ (/io) "~ > — H an d Cg (yUo) - ^ — fj which is equivalent to including only charm contribu- 
tions to the matching conditions for the corresponding operators. Analogously, only top loops 
contribute to X t . The last term in eq. (|2.2| ) now appears in X c , because it is entirely due to 
effects of charm loops in the RGE evolution. The splitting of charm and top is performed at 
the level of SM Feynman diagrams, and the effective theory is nothing but a technical tool for 
resumming large QCD logarithms in gluonic corrections to those diagrams. 

i - V 

0.2 0.4 . _ - - .0.6- ----- -0.8 . - 1 

-0.1 - - - 

-0-2 - . 

-0.3 - 

-0.4 - , ' 

-0.5 - 

-0.6 - ~ - 

Figure 2: X c as a function of rj (solid line), and its three components in eq. (|2~5| ) (dashed lines). 

X c is a function of rj that varies very slowly in the physically interesting region 0.4 < r) < 1. 
This is illustrated in fig. |4], where the three components of X c in eq. ( |2.5| ) are plotted as 
well. The second component is numerically small, while there is a strong cancellation of the 
//-dependence between the first and the third component. However, these components are not 
separately physical in any conceivable limit, so the cancellation cannot be considered accidental. 

Since X c is practically scale-independent, X t must be the source of the factor of ~3 en- 
hancement of BR 7 by QCD logarithms. This is indeed the case, because all the powers of r\ in 
eq. ( f2.6|) are positive and quite large. When r\ changes from unity to 0.566 (which corresponds 
to hq = My/ and fib = 5 GeV), then X t decreases from 0.450 to 0.325. At the same time, X c 
changes by only 0.008 (from — 1| ~ —0.639 to —0.631). Consequently, |C7°^ eff (/x)| 2 increases 
from 0.036 to 0.094, i.e. the branching ratio gets enhanced by a factor of 2.6. 

It is easy to identify the reason for the strong ^-dependence of X t . It is the large anomalous 
dimension of mb(fj,) that stands in front of the operator P-j ( |2.1| ). The anomalous dimension 
7 m is responsible for || out of || in the power of rj that multiplies the (numerically dominant) 
function A^x) in the expression for X t . Thus, the logarithmic QCD effects in b — > 57 can be 
approximately taken into account by simply keeping m;, renormalized at /^o ~ {m t or My/) in 
the top contribution to the decay amplitude. In section [|, we shall see that those features carry 
on at the NLO, lending support to the idea that they are not due to a numerical coincidence, 
and remain valid at the NNLO, too. 
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It would be interesting to understand the physics behind the different optimal normaliza- 
tions of m b in the charm and top loops. In this respect, the following observations may be 
helpful. Off-shell b — > S7 and b — > s gluon diagrams mediated by charm loops give rise only to 
dimension-six operators P aFD = -^s L {a^F u ^ , 0}b L and P aGD = -^s L {a f , u G aufl T a , 0}b L , 
respectively. The basis of physical operators in the off-shell effective theory for the charm sector 
can be chosen to be {P\, Pq, P<jfd, PcGd}- The LO Wilson coefficient of P fd is given by 
eq. Q2.5|) , i.e. its RGE evolution is very slow. When the evolution is terminated at //&, and the 
equations of motion (EOM) are used afterwards, P g fd reduces to the operator P 7 that contains 
an overall factor of m b . Thus, the b mass naturally associated with charm loops is a low energy, 
low virtuality mass. 




(a) (b) (c) 



Figure 3: Chirality flows in diagrams contributing to C\ (po)- The photon couples to any of the 
internal lines. 

In the top loops, on the other hand, mj, appears in two ways: either through the same 
mechanism as in the charm loops (fig. |3]a,b), or via the bottom Yukawa coupling y b of a right- 
handed b with a left-handed top quark (fig. |3|c). In the first case, the RGE evolution of C* FD (//) 
and C* GZ) (/i) is much faster than in the charm sector. Such a fast running can be compensated 
to a large extent by using rrib(fio) in the EOM at the low-energy scale. In the second case, it is 
not necessary to use the EOM to project on P 7 . The natural normalization scale of the Yukawa 
coupling yb{n) is given by the heavy masses circulating in the loop. The b mass associated 
with the Yukawa contributions is therefore a high-energy mass, m&(/io) with /x ~ iMw or m t)- 
This reasoning leads us to the conclusion that the appropriate normalization of the total top 
contribution is in terms of m&(/Xo). 

A careful reader might worry about how this picture can be reconciled with GIM cancellations 
that take place in the case of degenerate quarks. Some light on this issue can be shed by 
considering a hypothetical situation when m b <C m c <C m t <C M w . In such a case, the effective 
6S7 interaction would be relatively well approximated by 

AC ( 23 1 

£ int « -^V*V tb {-KW - m b (m c )] + 0(m c 2 t /M^)| s L a^b R F^. (2.8) 

It is not an exact formula even at LO, because we have neglected the very slow RGE running 
from My/ to m c t and the small anomalous dimensions of P 7} %/m b (n). However, eq. ( J2.8|j works 
numerically quite well. The QCD enhancement of BR 7 can be much larger here (than in the real 



7 



world), and its main reason is clearly seen. The number || arises from the matching condition 
with an effectively massless internal quark. In the absence of QCD, it drops out owing to GIM 
cancellation. However, no cancellation takes place any longer when the QCD evolution of m b 
is taken into account. 



3 Next-to-leading order formulae 

We have seen in the previous section that most of the leading QCD logarithms can be taken 
into account by keeping the top and charm contributions split, and by renormalizing m&(//) 
in P 7 Q2.1p at jiQ ~ {m t or M w ) in the top contribution to the decay amplitude. Motivated 
by this observation, we shall now rewrite the known NLO expressions for B — > X s ^f in an 
analogous manner. As we shall see, this simple operation not only allows us to reproduce the 
logarithmic QCD enhancement in a satisfactory manner, but also the NLO corrections become 
significantly smaller than in the traditional approach. Moreover, the residual renormalization- 
scale-dependence diminishes, without any accidental cancellations involved. In other words, 
the behaviour of QCD perturbation series improves. 

Assuming that the dominant NNLO QCD effects have the same origin as the dominant LO 
and NLO ones, we shall use all the currently known perturbative information to determine 
the ratio of mfc(/i ) to the low-energy 6-quark mass that normalizes the semileptonic decay 
rate. This low-energy mass will be chosen in a manner that ensures good convergence of 
Q CD-perturbation series for the semileptonic decay. 

Our input here are the standard NLO QCD formulae for B — > X s 7 collected in ref. |T2|| . 



and the separate charm-sector and top-sector matching conditions for the relevant operators 



presented in section 2 of ref. ||23|| . The evaluation of such separate matching conditions was 
described in great detail in section 5 of the latter paper. We have also checked it independently 
with the method of ref. [IJ. Apart from the perturbative QCD effects, we shall include the 
electroweak and the available non-perturbative corrections. 

The B — > X s 7 branching ratio with an energy cut-off E in the 5-meson rest frame can be 
expressed as follows: 

v t * s v tb 



BR[B -> A s7 ]?>C d *' * = BR[-B X c ev\ 



cxp 



V, 



cb 



" m [P(E ) + N(E )}, (3.1) 



7T 



c 



where a em = a™ she11 |T7|] and P(E ) is given by the perturbative ratio 



rj> -> x s ^] Ei>Eo 

IWK 6 | 2 r[6^ X u ev\ 



v; s v tb 



(,n "' n P(E ). (3.2) 



7T 



N(Eq) denotes the non-perturbative correction.^ Contrary to the standard approach, we have 
5 This means that P(Eq) gets replaced by P(Eo) + N(Eq) when b is replaced by B in eq. 



chosen the charmless semileptonic rate (corrected for the appropriate CKM angles) to be the 
normalization factor in eq. (|3.2|) . This modification is offset by the factor C in eq. ( [3.11) : 



C 



V, 



ub 



F[B -> X c eV] 



(3.3) 



V cb r[B^X u eu) 

This observable can either be measured or calculated. Our normalization to the charmless 
semileptonic rate in the l.h.s. of eq. ( |3.2| ) is motivated by the need for separating the problem 
of m c determination from the problem of convergence of perturbation series in b — > Xgj. The 
factor C can be called "the non-perturbative semileptonic phase-space factor" . 

The superscript "subtracted if) and ip'" on the l.h.s. of eq. (|3.1| ) means that the processes 

B -> if)^xW 

\\T (2) 7 (3.4) 

are treated as background, and should be subtracted on the experimental side. This background 
is negligible (below 1%) for a cut-off of 2.1 GeV on the photon energy in the S-meson rest 
frame, but may have a ~ 5% effect when the cut-off is lowered to 1.8 GeV.[] 
The perturbative quantity P(Eq) can be written in the following form: 



P(E ) 



K c + 



1 H In r(fi )K t + e e 

7T mf 



(3.5) 



where K t contains the top contributions to the b —>■ s'j amplitude. K c contains the remaining 
contributions, among which the charm loops are by far dominant. The electroweak correction 
to the b — > S7 amplitude is denoted by £ ew . The ratio 



m 



MS 



1 f 



(3.6) 



appears in eq. fl3.5p because we keep m b renormalized at /i in the top contribution to the 
operator P 7 ( |2.1| ), while all the kinematical factors of m& are expressed in terms of the bottom 
"IS mass". The IS mass is defined as half of the perturbative contribution to the T mass. As 



argued in ref. |2T|, expressing the kinematical factors of in inclusive 5-meson decay rates in 
terms of the IS mass improves the behaviour of QCD perturbation series with respect to what 
would be obtained using m^ s (m 6 ) or ml° le . 

The bremsstrahlung function B(Eq) contains the effects of b — ► sjg and b — > sjqq 
(q = u, d, s) transitions. It is the only ^-dependent part in P(E ). It is given in appendix E. 
Its influence on the b — ► X s 7 branching ratio is less than 4% when 1 GeV < Eq < 2 GeV. 
Therefore, we do not split the top and charm contributions to this function. It would not 
improve the overall accuracy at all, but only make the formulae unnecessarily complicated. 

6 The effect would become more than 100% for E - 0.2 GeV. 
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16 
23 
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23 


12 
23 
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0.4086 


r\ /loon 

—0.4230 


—0.8994 
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1.41U f 


— U.oooU 


n a oo^ 
— U.4/OD 


— U.U i 14 


— U.b4y4 


n noon 
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— U.UUo 1 


di. 


— 17 6507 


11 3460 


3 5762 


—2.2672 


3 9267 


1 1366 


— D 5445 


1653 




9.2746 


-6.9366 


-0.8740 


0.4218 


-2.7231 


0.4083 


0.1465 


0.0205 










1 


1 






















1 


1 
















0.4702 





-0.4938 


-0.4938 


-0.8120 


0.0776 


-0.0507 


0.0186 


e k 


5.2620 


-3.8412 








-1.9043 


-0.1008 


0.1216 


0.0183 



Table 1: "Magic numbers" for K c and K t 



Our result for K r reads 



K c = J2v ak \d k + 
k=i I 



4tt 



2f3 Q a k d k f In — + 77 In 



+4 + d v k 7] + d a k a{z) + d\b{z) + cf/ivr 



} + ^S ( ^ + ^ 4)[aW + 6W] ' (3 - 7) 



where (3q = y and z = (m c /mi,) 2 . The functions 0(2;) and are given in appendix D. 

The "magic numbers" d k , d k , ... can be found in table [1]. In their evaluation, the unknown 
NLO matrix elements r 3 , ...,r 6 [r2~| have been set to zero. The resulting error will be absorbed 
in the NNLO uncertainty below. This is mandatory^] because C^ 1 ' / (l^b) < for 
k = 3, 6. 

The NLO expression for K t is as follows 



K t 



1 - -a s (m b ) + 
y 



Q! s (A*o) 



d 



71 



In — 4x— — 

m t ox 



1 _L .f / \ 4/J^ _2_ \ j / , 

--7723 + - - 7723 j F { (x) 



47T 



I fc=l 



4-^23 



7411 



1 . t , , / 12523 

—riAUx) + 

2 / n ; 1 3174 47gl 



n 



+-tiFHx) + , 
3 ' 1V ; V 4761 



50092 1110842 
+ 



-7T 

9 
16 



77 H tt" H I In h 77 In — 

357075 '27 9 \ fi b ' m t 



\ (ln^ + ,,ln«)U(*) 
3 \ J J 

F*(x) 



32 



+7723 



4 . /2745458 38890 4 , ., 16 / m b , /i \\ , 

— 77F|(x) + 77 - -Trfvr + 7) In — + 77 In — F^x) 

3 ' 1V 7 I 357075 14283 ' 9 V ' 9 1 ' ■ 11 " 



fib 



m t 



(3i 



We have checked that those parts of r%, ...,tq that are not due to closed 6-quark loops have only a 0.5% 



effect on BR^ 
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The functions Aq and Fq of x = (m t ([io) / Mw) 2 have already been given in eq. (|2.7|). The 
remaining functions read 

zpt{~\ — -9a 2 +16a:-4 i„ , 7z 3 +21x 2 -42:z-4 
- 6(x-l) 4 UlX-f- 36(x-l) 3 ' 

/it /_\ _ 32x 4 +244x 3 -160a; 2 +16a; r • /i 1^\ , -774x 4 -2826x' 3 +1994x 2 -130:r+8 i_ _ 
A iy X > ~ L ^\ i ~x) + 81(l-a0« lnX 

-94x 4 -18665x 3 +20682a: 2 -9113x+2006 



+ 



243(l-a;) 



4 • 



T?t/„\ _ 4x 4 -A0x 3 -41x 2 -x j ■ {-, l\ | -144:r 4 +3177:r 3 +3661x- 2 +250:r-32 1v , „ 
- SO^) 3 2 V l ~x) + 108(1-1)5 

I -247x 4 +11890x 3 +31779a- 2 -2966^+1016 
^ 648(l-x) 4 • 

The 0(a 2 ) term is included in K t along the lines of the T expansion ^TJ. Its effect is at the 
level of 1% only, and it is offset by an analogous term in r ( |3.6|) (see appendix B). 
The electroweak correction e ew in eq. ( |3.5| ) consists of three terms 

e ew = 6™Ci 0)cS ^ b ) + ^^Cr (0)eS (» b ) ~ (K^ + rKi 0) ) In (3.9) 

a s (//bj 7T v 'lib 

The first term stands for the dominant non-logarithmic electroweak effects calculated in ref. [[0| 
(see eq. (11) of that paper). The remaining two terms originate from the logarithmically en- 
hanced electromagnetic corrections to the weak radiative (eq. (12) of ref. |18|]) and semileptonic 
(eq. (1) of ref. P^f| ) decays, respectively. In the last term, Kf^ and are just the LO 
contributions to K c and K t . 

The non-perturbative correction N(Eq) in eq. ( |3.1| ) is partly known 

N(E ) = -1 (X CO) + rK f) + rri) ^ + - (3.10) 

The calculable A 2 /m 2 correction is taken into account above, while the calculable A 2 jm\ 
corrections [§] have cancelled out because of our normalization to the charmless semileptonic 
rate. The dots stand for higher-order terms in the heavy-quark expansion |7|,[?]J, as well as for 
the non-perturbative effects due to higher (than if) and if)') intermediate cc states, to light-quark 
loops and to motion of the 6-quark inside the B meson. In our numerical results (calculated with 
E = 1.6 GeV), we shall set those effects to zero without including any additional uncertainty. 



11 



4 Numerical analysis 



In the present section, we shall test our formulae numerically. The experimental inputs and 
the ranges for the renormalization scales are collected in appendix A. Appendices B, C and D 



are devoted to the determination of the mass ratio r ( |3.6| ), the phase-space factor C (|3.3| ) and 
the z-dependent terms in K c (|3.7|) , respectively. The final results obtained there read 



r(p = m t ) = 0.578 ± 0.002^ ± (parametric errors) (4.1) 

C = 0.575 (1 ± 0.01 ±0.02 ±0.02) (4.2) 

a(z) = (0.97 ±0.25) + i(1.01 ± 0.15) (4.3) 

b( z ) = (-0.04 ±0.01) + i(0.09±0.02). (4.4) 

In our NLO computation, the imaginary parts of K c and K t are irrelevant, because all the 
0(a 2 ) terms on the r.h.s. of eq. (|3.5|) are set to zero after the square is taken. If the imaginary 
parts of K c and K t were not set to zero, they would affect P(E ) by only 0.5%. 

The ratio of CKM angles standing in front of P(E ) can be expressed in terms of the 
Wolfenstein parameters as follows 



v t *v tb 2 



V, 



cb 



1 + A 2 (2p - 1) + A V + f- A 2 ) ± 0(X e ) 

0.971 + 0.10(p - 0.224) = 0.971 ± 0.004, (4.5) 



where we have used A = 0.2237, A = 0.819, p = 0.224 ± 0.038 and fj = 0.317 @. The only 
relevant source of uncertainty is the error in p. Even if this error were enlarged by a factor of 
2, the influence of | V t *Vtb/V c b\ 2 on the overall uncertainty in B — > Xg'j would remain negligible. 
The central value of eq. (|4.5|) is also consistent with the analysis of ref. [27|. 



The electroweak and non-perturbative corrections from eqs. ( |3.9| ) and ( |3.10| ) take the follow- 
ing values^ 

e cw w 0.0035 ±0.0012 + 0.0028 = 0.0075 (4.6) 

N(E ) = 0.0036 ±0.0006. (4.7) 

Their effects on BR 7 are —3.8% and ±2.5%, respectively. In the electroweak correction, 
^Higgs = 115 GeV has been used. The sensitivity of our final result (|4.14|) to Mmggs is very 
weak (only 0.3% when M ffiggs changes from 115 GeV to 200 GeV). In N(E ), A 2 = 0.12 GeV 2 



and m c = m c (m c ) = (1.25 ± 0.10) GeV |29| have been used. The indicated uncertainty is due 
to m c only. 



8 The electroweak correction remains practically unchanged (e cw = 0.0071) when the very recent results of 
ref. [E8fl are included. 
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In the "naive" approach, when only the one-loop electroweak diagrams are calculated, one 
finds 

23 1 

K n^ vc = w _ .639 anc i jq naiv0 = — Al(x) « 0.450. (4.8) 

36 2 

From those numbers, we obtain 

K™ ive + mb ^ K™ ivc « -0.379, (4.9) 

which implies BR" aive « 3.53 x 10~ 4 (for N(E ) = 0). If the factor m b (m t ) / 'm\ s were not 
included, our "naive" prediction would be 4 times lower .0 



At LO, i.e. when the 0(a s ) terms in eqs. ( |3.7|) and fl3.8|) are neglected, we find 



tff = -0.631SS and = 0.434+^1 (4.10) 



which implies that BR^ = 3.56jlo;o7 x 10 4 . The indicated errors correspond to the variation 
of the low-energy scale as described in appendix A. 
At NLO, we find for E = 1.6 GeV 

K c = -0.61llS + i (-0.032^), (4.11) 
K t = (0.397 ± 0.003) +i (0.011 ±0.002), (4.12) 
S(Eb) = 3.11^ x 10~ 3 . (4.13) 

From the above three results, we obtain BR[_B — > X s ^] Ei> i q GeV = 3.601o;q 4 x 10~ 4 . As before, 
the quoted errors are due to /x& only. 

We can see that the //&- dependence of Re(i^ c ) and R,e(K t ) is very weak, both at LO and 
at NLO. Such a weak scale-dependence is somewhat surprising at LO, but at NLO it is not. 
Moreover, at NLO, the weak scale-dependence is not a result of any accidental cancellation 
among strongly dependent terms. Whatever cancellations occur, they occur inside the charm 
contribution K c , and thus cannot be considered accidental. 

Since the dependence of the NLO branching ratio is very weak (below ±1.4%), and the 
dependence on no is not stronger, we can estimate the NNLO corrections by simply saying that 
they are of order (a s (mb)/7r) 2 w 0.5%, up to an unknown factor of order unity. Thus, it seems 
safe to assume a value of ±4% for the theoretical error which is due to neglecting r%, 
at NLO and to the NNLO effects. This is almost twice the combined scale-dependence of the 
result obtained by scanning /x& and /io between one half and twice their central values (2.2%). 
We consider the error related to the value of m c /mb in K c separately. The latter uncertainty 
amounts to ±5.5% alone, when our estimate from appendix D is used. 



If m b (M w )/m b (5 GeV) « 0.713 were used instead of m b (m t )/m^ s » 0.578 in P{E ) Q3.5|), we would 



obtain BR™™ w 2.49 x 10~ 4 . 
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When all the sources of uncertainties are included, we find 



BR[B -> X^l^tev * = 3.60 x 10" 4 (1 ± 0.06 (mc/mt in Ko) ± 0.04 (other nnlo) 

±0.01 (pcrt C )±0.02 Al ±0.02 A 

±0.02 (M z ) ± 0.02BR(semilept) cxp ± 0.01 mt ) 

= (3.60 ± 0.30) x 10" 4 . (4.14) 

The errors in the second line above originate from the semileptonic phase-space factor in 
eq. ( C.13 ). The last of them is due to the error in m\ s . Other effects caused by the uncertainty 
in m\ s are negligible (below 0.2%). 

In the last line of eq. ( [4.14| ), all the uncertainties have been added in squares. Thus, the final 
error has only an illustrative character, because many of the partial ones have no statistical 
interpretation. 

We have chosen Eq = 1.6 GeV instead of the commonly used Eq = ~ 0.23 GeV (i.e. 

5 = 1 — 2E /mb — 0.9), because the low-energy photons in B — > X s 7 are not experimentally 
accessible. For high values of the cut-off (around and above 2 GeV), non-perturbative uncer- 
tainties are much larger because the S-meson shape functions are unknown. Models for the 
shape functions considered in refs. |TT],[J0] suggest that E = 1.6 GeV is already low enough to 
make those uncertainties negligible. In our opinion, the way to resolve this issue is to measure 
the photon spectrum without any theoretical input, with the photon energy cut-off as low as 
possible (e.g. 2.0 or 1.9 GeV in the 5-meson rest frame). Next, a simple extrapolation to the 
theoretically known spectrum in the vicinity of 1.6-1.7 GeV should be made. For this purpose, 
the following approximate expression for the integrated branching ratio as a function of Eq 
might be useful: 

BR(Eq) = (3.524 + 0.298Eo - 0.157£ 2 ) x 10~ 4 , (4.15) 

where Eq is expressed in GeV. The above expression has been obtained as a fit to our results 
in the region 1 GeV < E < 2 GeV. The accuracy of this fit is ±0.3%. Outside the considered 
region of Eq, our formulae from section § are not expected to work well. On the low-energy 
side, it is due to the fact that we have not included the fragmentation functions discussed in 
ref . || . As far as the high-energy side is concerned, the growth of uncertainties with Eq can be 
roughly estimated by using fig. 3 of ref. |10| and fig. 3 of ref. | 30fl . However, translating those 
plots into quantitative estimates is rather difficult, because of the model-dependence involved. 
A further study of this issue is crucial for a precise comparison of theory and experiment in 
B -> X sl . 

In view of the fact that many of the published results have been calculated for S = 0.9 (i.e. 
Eq = ^mb), it is interesting to check what our formulae give in such a case. We find 

B (±m b ) = 8.2± 3 2 l x 10~ 3 . (4.16) 
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Together with eqs. (|4.11|) and ( |4.12 ), this gives 
BR[B -> X sl ] Ei>mb/w = 3.73 x 10" 4 . 



(4.17) 



The total relative error in this case is roughly comparable to the one given in eq. ( 4.14j) . If 
we used m c /mb = 0.29 instead of 0.22 in K c and B(Eq), we would find 3.35 x 10~ 4 for the 
branching ratio. The latter result is very close to the ones obtained in many previous analyses 
(see e.g. |H|[I9|]). Thus, the replacement of m^ ole /ml° le by m^ IS (^)/m^ 5 in (P 2 ) is the main 
reason why our result is significantly higher than the previously published ones. 

Our NLO formulae differ from the ones used previously by residual 0(a 2 s ) terms. Quite 
unexpectedly, those terms tend to cancel among themselves, even though they are not individ- 
ually very small. For instance, if we used eq. ( P-l|) instead of (|B.4| ) for the determination of the 
ratio r, our final result for the branching ratio would be 7% lower than the one in eq. ( f4.14p . 



5 Constraints on new physics 

In the absence of new light degrees of freedom, physics beyond the SM manifests itself through 

(i) new contributions to coefficients of the operators involved in the SM calculation and 

(ii) the appearance of operators absent in the SM, such as operators with different chiral- 
ity [^TJ. Examples of contributions beyond the SM include diagrams like the ones in fig. 1, with 
charged Higgs exchange or with chargino-squark loops. Therefore, new physics contributions 
are at the same level as the SM ones, and similar theoretical accuracy is necessary. 

In all the new physics scenarios which do not involve new operators (and in which 
C£ ew (/i ) =0 for k = 1,2,3,5,6), it is straightforward to incorporate extra contributions 
to the NLO formulae of section ||[ These contributions effectively modify eqs. flpDQ and (E.9) 
through the replacements 

A){x) - 4(*)-2Cf )new (/i ), F*(x) - i^(*)-2Cf new M, (j = 0,l) 
El(x) - ^(aO+Cl 1)ncw (/i ), (5.1) 

where C^ new are the LO (j = 0) and NLO (j = 1) new physics contributions to the Wilson 
coefficients of the operators Pi 

CrVo) = )n ° w (/io) + ^ (1)nCW (^o). (5.2) 

47T 

Beyond the SM, the necessary NLO coefficients are known at present in the general Two- 
Higgs-Doublet-Model (2HDM) in some specific SUSY scenarios |, |, || |35[ , and 
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Except for the terms proportional to ln(/Jo/mt) in eq. (|3.8| ). Such terms should be left unchanged, because 



the appropriate logarithms of fio are already contained in C7g ew (/io). 



15 



in the left-right symmetric models | 53fl . The coefficients in ( |5.1| ) are evaluated at the scale 
Ho ~ M w ,m t . 

Whenever the new particles have masses of order A ^> My/, they should be decoupled at 
/i ~ A. The coefficients C^ new should be then evolved down to the scale /i , in order to resum 
all large logarithms correctly. Examples of this procedure can be found in refs. 



36] 



Notice that the coefficients in eq. (|5.1|) are not the effective ones [22], in terms of which the 
matching conditions are often given. The relation between the two sets is very simple: 



(5.3) 



where y^j = (0, 0, — | 



3' 
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y), and all the other vanish. 
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Figure 4: Direct and indirect lower bounds on M H + from different processes in type II 2HDM as a 
function of tan/5. The B — > X s ^ bound is the one in eq. ( |5.5| ) below. 

To illustrate the incorporation of new physics in our calculation of B — > X s 7, we consider 
the case of the 2HDM, which also exemplifies well the importance of this decay mode for new 
physics. We shall update the lower bound on the charged Higgs mass in the type II model. 

In the 2HDMs with vanishing tree level FCNCs, the only additional contribution with respect 
to the SM comes from the charged Higgs boson-top loops. It depends on the mass of the charged 
Higgs boson, Mh , and on the ratio of the v.e.v's of the two Higgs doublets, tan (3. Models of type 
I and II differ by the way fermions couple to the Higgs doublets. In the type II model (realized 
in the MSSM), the charged Higgs loops always enhance BR 7 , while the decoupling occurs 
slowly. Therefore, this decay mode provides strong lower bounds on Mh, whose dependence 
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on tan/3 saturates for tan/5 ~ 5. Previous calculations led to Mh ^ 250 GeV, independently 
of tan/3 Hlj[32"l, 37j. This bound is much stronger than the one from direct searches at LEP2 
(Mh > 78.5 GeV [@), and than the indirect lower limits from a number of other processes 
(see fig. [|) |5^,^D|. In model I, B — > X s ^ is less restrictive than other processes, because the 
charged Higgs loops tend to suppress the branching ratio and decouple for large tan/3. The 
most important constraint in that case comes from Rt> [|TJ. 

The LO and NLO Wilson coefficients in the 2HDM are given in eqs. (52)-(64) of the first 
paper in fL]. Adopting the same notation, with y = (m t (fi ) / M H ) 2 and Hw = Ho, we perform 
the following replacements in K t 



El{x) -> E t (x)+E H (y), 



A\(x) 
F{{x) 



A\{x)-2 
Fl{x) - 2 



G?(y) + A?(y) In 
Gf (y) + Af (y) In 



f4 



M H 



F t (x)-2^ 0)eff ( /Uo ) 



(5.4) 



As mentioned above, the In /i terms in eq. ( |3.8| ) should be left unchanged. 
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Figure 5: The 99% CL bound on the 2HDM-II charged Higgs mass from B — > X s ^ as a function of 
the world average and of its error. The contour lines represent values which lead to the same M H + 



bound. The experimental world averages evaluated with use of the preliminary [fTjJ and published 
CLEO results are indicated for reference. 
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We are now in a position to calculate the branching ratio for B — > X s ^ in the 2HDM, and to 
find the lower bound on M# in model II. Since the CLEO and BELLE results (obtained with 
Ky > 2.0 GeV and 2.1 GeV, respectively, in the T frame) have been extrapolated to the "total" 
rate, we compare the experimental result with our prediction for > ^rrib JTUJ. With respect 
to previous analyses, the SM prediction is now higher — see eq. Q4.17 ). Moreover, the charged 



Higgs contribution in model II cannot help reducing the value of BR 7 . Thus, our bound is 
going to be stronger than in the past [|l],|32],[37]]. In addition, one should take into account that 
for a heavy charged Higgs, the dependence of BR 7 on Mh becomes very mild, signalling the 
decoupling. Consequently, a small change in BR 7 affects the M H bound in a significant way. 

Therefore, we adopt a very conservative approach. We employ the experimental world 
average (3.23 ± 0.42) x 10~ 4 and scan over the main theoretical errors in our calculation - 
first line of eq. ( f4.14|) — notably the one related to m c /mi, in (P2). Whenever the and 
/^o-dependences of BRy in the 2HDM are larger than our 4% NNLO uncertainty in eq. ( [4.14 ), 



we expand this error accordingly^] All the remaining parametric and experimental errors are 
combined in quadrature. The absolute lower limit is 

M H > 350 GeV (5.5) 

for any value of tan/3. As far as the gaussian errors are concerned, this is a 99% CL bound. 
The 95% CL bound is 500 GeV. In view of future changes in the experimental situation, we 
display in fig. |5| the 99% CL bounds as functions of the central value and of the error of the 
world average. The dependence of the Mh lower limit on tan (3 is shown in fig. |j. 

The bound is clearly very sensitive to the way various errors are combined. For instance, if 
we combined all the errors in quadrature, the 99% CL absolute lower limit on Mh would be 
500 GeV. On the other hand, the way we combine the errors emphasizes the large uncertainty 
coming from m c (ii)/m\ in {P%}- In effect, the lower limits that we have quoted correspond 
to m c {m c ) / m\ s = 0.26, i.e. the highest value of m c {ji)/m\ s that is compatible with our 
analysis in appendix D. Had we employed m c /rrib = 0.29 ± 0.02 as in previous analyses and 
treated this error as in the derivation of eq. (|5.5[) , the 99% (95%) CL bound would have been 
M H > 280 (380) GeV. 

In deriving the limits on Mh, we have not assumed that the 2HDM is valid for sure, i.e. we 
have not subtracted xl CEtfit from the \ 2 used for the test. Consequently, the statement that 
M H > M at x% CL means that the 2HDM is excluded at x% CL unless M H > M . A different 
procedure^ consists in assuming the validity of the 2HDM and in testing x 2 ~ x\ ■ In this 
case, adopting the same treatment of errors as in eq. ( |5.5| ), we obtain M H > 345 (490) GeV at 
99% (95%) CL. 



11 Here, the scales are varied between 0.4 and 2.5 times their central values. 

12 We thank A. Strumia for bringing this point to our attention. 
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6 Conclusions 



The reanalysis of B — > performed in the present article contains several major improve- 
ments with respect to the previous ones. In particular, the semileptonic phase-space factor 
is expressed in terms of an observable for which the NNLO results are known. Moreover, the 
RGE evolution of in the top-quark contribution to the b — > 37 amplitude is identified 

as the main reason for the huge enhancement of the branching ratio by QCD logarithms. In 
consequence, better control over the behaviour of QCD perturbation series is achieved. Unfor- 
tunately, the desired reduction of theoretical uncertainties significantly below the current exper- 
imental ones is found to be impossible at the NLO level, because of the strong m c -dependence 
of certain two-loop diagrams with charm-quark loops. 

In order to compare our final prediction (|4.14j ) with the weighted average ( |1 . 1|) of the exper- 
imental results, we have either to evaluate the theoretical errors in eq. ( f4.17| ) or rather to "take 
back" the extrapolation to low photon energies on the experimental side. Choosing the latter 
possibility, we find that the experimental result corresponding to our eq. ( [4.14 ) is 

BR[B - X a7 ]^>i.6 Gev = f^( 3 - 23 ± °- 42 ) x 1(r " = ( 3 - 12 ± °' 41 ) x 10 " 4 > 

where the ratio of eqs. ( f4.14p and ( |4.17 ) has been used as the rescaling factor. The difference 
between theory and experiment can now be written as 

BRf - BR^ xp = (3.60 ± 0.30) x 10" 4 - (3.12 ± 0.41) x 10~ 4 = (0.48 ± 0.50) x 10" 4 , (6.2) 

which is compatible with zero. Based on the present experimental information, we have shown 
that the Two-Higgs-Doublet-Model II with a charged Higgs boson lighter than 350 GeV is 
strongly disfavoured. 
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Appendix A 



29[, we take 



In this appendix, our numerical input parameters are collected. From ref. 
M z = 91.1882 GeV, M w = 80.419 GeV, a em = 1/137.036, a s (M z ) = 0.1185 ± 0.0020, 
rof(mf) = (1.25 ± 0.10) GeV and BR[B -> X c ez/] cxp = 0.1045 ± 0.0021. 
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For the 6-quark mass in the so-called lS-scheme, we use m\ s = (4.69 ± 0.03) GeV ||4"2| . As 
far as the top-quark mass is concerned, mf IS (mf IS ) = (165 ±5) GeV is used, which corresponds 
to m^ ole = 174.3 ±5.1 GeV 



In several places, definitions of m t and m h have been left unspecified, because it would 
become relevant only at higher orders of perturbation theory. In such places, m^ [S (m^ [S ) and 



m\ s are used. 



When the /^-dependence is being tested, each scale is made to vary from half to twice its 
central value. The central value of ^ is chosen to be m b . For the matching scale /i , we take 
the central value of Mw in K c and in the quantities for which no flavour splitting is performed. 
In K t and in the mass ratio r(fi Q ), the central value of fiQ is set to m t . 

The final result (|4. 14 ) for the branching ratio ( |3.1|) has been found by using the above inputs 
as well as the quantities specified in eqs. (pD, flOD , ( gTTD , (|[|), (|UT3|) and ( [TXBl) . 



Appendix B 

Here, we determine the mass ratio r(/z ) = m b (fi ) MS / ml s . The NLO result reads 



>(>o) 




7462 a a (fi ) 15926 



1587 



+ -a s (m b ) 2 



0.611. 



fiO=m t 



(B.1) 



1587 a s (mb) 

The 0(a1) term is included along the lines of the T expansion Its effect on r is at the 
level of 1% only, and it is offset by an analogous term in K t Q3.8Q . 

Since m;,(/io) could, in principle, be determined from a high-energy observable (e.g. the 
Higgs width), we are allowed to calculate the ratio r(/i ) using as many orders of perturbation 
theory as are currently known, i.e. up to 0(e 3 ) in the T expansion. 



From eq. (168) of ref. | I2| , we find 



MS, 



m MS) 



15 



0.889 



'as(Mz) 
0.1185 



-0.26 



ni 



1S\ 0-04 



4.69 



m c {m c 
1.25 



-0.003 



4.69 



0.006 



(B.2) 



All the masses and renormalization scales in eqs. ( [B.2|) -( p74] ) are expressed in GeV. The de- 



pendence on /ib in eq. (p.2| ) originates only from the fact that a s (nb) rather than a s (rrib) has 
been used in evaluating this relation. Varying /i^ by a factor of 2 around 4.69 GeV gives us an 
estimate of ±0.4% for the neglected higher-order effects. 
The four-loop RGE for the quark mass imply 



m 



MS, 



m t 



° v Jl w 0.650 



m MS( m MS) 



'as(M z ) 
0.1185 



-0.72 



m 



MS, 



m 



0.19 



4.17 



U65 



-0.08 



(B.3) 
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where the higher-order uncertainty is negligible (below 0.2%). Combining eqs. 
we find 

-1.0 / i <?\ 0.23 / / \ \ -0.003 



and Q , 



r{m t ) = 



m 



is 



0.578 



a s (M z ) 
0.1185 



< l e\ 0.23 

m b \ 
469 J 



m c {m c ) 
1.25 



U65/ 



0.08 



4.69 



0.006 



(B 



The above approximate formula works with better than 0.2% accuracy (compared with the 
complete one0 ) when the input parameters vary within their errors (see appendix A) . 

It is important to note that the central value in eq. ( p. 4j ) is 5.4% less than r NLO in 
eq. ( |B.1| ). The main reason of this suppression is the —3.2% higher-than-NLO correction to 
l )/m\ s (see eq. (99) of ref. f4"2"f). Switching from two-loop to four-loop RGE causes a 



m 



MS^ m MS 
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reduction of r by only 0.7%. Additional effects originate from the fact that the scales /i = m i 
and = m b (m b ) have been identified in eq. (|B.1| ). 

The ratio r(/io) depends only logarithmically on the actual value of m\ s . This is reflected in 
eq. ( p.4|) where the power of m^ 5 /4.69 is equal to 0.23 only. Consequently, our results are not 
very sensitive to the numerical input we use for m\ s . They would not change much if instead 
of m\ s we used another definition of the "kinetic" mass of the b-quark, e.g. the ones defined 



in refs. |43|,^4|, or the MS mass at appropriately chosen fi. However, in such a case, eq. 
would need to be modified accordingly. 



Appendix C 



The present appendix is devoted to a determination of the semileptonic phase-space factor: 
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(C.l) 



along the lines of the T expansion [^T 

We begin with t 
masses and a s = a 

T[B 



We begin with the NNLO expressions for both decay rates expressed in terms of pole quark 
(n,=5 \m b ) ^0.219: 
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Vcb\ g{z P 



l +e %(D (z)+e 2^(2) ( 
7T ' 7T Z 




6(1 - z n 



2m 2 b 



1927T 3 



V ub \ 2 



l + e%«+6 2 ^(^) + 

7T TX Z 



Ai 
2m 2 b 



9A 2 
1m\ 



(C.2) 
(C.3) 



where z n 



[m 



polc /m^ olc ) 2 and g(z) = 1 — 8z + 8z 3 — z A — 12z 2 lnz. The auxiliary parameter 



1 will be helpful below in denoting orders of the T expansion. 



13 



We thank A. Hoang for sending us the numerical data from which eq. (168) of ref. [[f2| was derived. 
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The perturbative corrections follows PS: 



P. 



«(*) 



0.09) 



2/i(z) 

-1.68/^ 4) + (1.4 ± 0.4) 



25 
If 



-7T 



11 



where fi 
h{z) = -(1- 

+(1- 
-(l + 16^ 2 + 2 4 )[6Li 2 (z 



-3.22/^ 4) + 5.53 + /„(*) 
and 



•12.6 ±0.4, 



-21.30 + f u (z), 



2 , /25 239 25 2 . , . 1( 
r ) ( —2 + — ) + z In * ( 20 



4 
17 



3 
64 



4 
17 



— - — z + — z 2 ) ln(l 
3 3 3 ' 



4 .2 , 17- 



90z - + —2' 
3 3 

4(l + 30^ 2 + ^ 4 )ln2 ln(l 



(C.4) 

(C.5) 
(C.6) 
(C.7) 

z 2 In 2 z (36 + z 2 ) 
-*) 



71 



32^ 3/2 (l + 



7T 



- 4Li 2 (v^) + 4Li 2 (-v^) - 2 In 



1 - 



1 + 

As far as the charm-mass correction f u (z) in eq. ( |CJ.7| ) is concerned, the only thing we know 
for sure is that /«(0) = 0. A rough estimate of / u (0.09) can be obtained by considering the 
relation between pole and MS bottom masses, for which charm-loop corrections are known at 



In z 



two [46 1 and three |^2| loops. Similarly to f u (z), they originate from charm loop insertions on 
gluon lines. For m c /mb = 0.3, we learn from ref. [32|] that two-loop charm mass effects lead 
to an extra term 0.43(a s /7r) 2 . As the charm loops are closely related to the BLM corrections, 
which amount to 1.56/?o(«s/ 7r ) 2 in the relation between masses, we can rescale them by a factor 
—3.22/1.56 to obtain / u (0.09) ~ —0.9. Of course, this is only a very rough estimate, and an 
actual calculation of f u {z) would be welcome. 

When the ratio ( |C.1| ) is calculated, the overall factors of (mf ole ) 5 cancel. The remaining 
explicit dependence on pole quark masses can be eliminated with the help of the following 
relation: 



rn 



pole 



rn 



pole 



m 



pole 



rn 



pole 



rn 



pole 



1 



m B - m D + 2=^ 



2m o 



A 3 



with Cf 



\m T (l 
9.460 GeV, m B = 



A 2 



{l + e^f^ [l + Y + 0{af)\ } 

; 5.314 GeV and 



(m-r/2) 2 

m T ~ y.4bU UeV, = |(3m B *± + m B ± + 3m B *o + m B o 
Wif) = |(3m D *± +m D ± + 3m D *o + m D o) « 1.973 GeV. In the denominator, A describes the non 
perturbative contribution responsible for the difference between |mx and m\ s . From ref. f4"2" 
one finds 14 
A 2 



(mr/2) 



(0.04 ± 0.03) GeV, 



(C.9) 
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Here, we neglect the 0(a i s ) difference between ai n/ 4 \rrib) and ai"' '''(nib) in the MS scheme. 



>/=5), 
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203 25 . , _ , 
Y = —-—\n(a s C F )+y 



3m c 
a s mr 



with 



2 a 7T 4 n 27T □ 

y(a) = -ln — I — a a H a 

yw 3 2 2 3 3 



4 -2a 2 



arctan • 



1 



(CIO) 



(C.11) 



Below, m c (m c ) « 



3Va2^T Va+T 

1.25 GeV [2J| will be used in this function. The influence of y(a) on C is less 
than 0.3%, so it does not matter what definition of m c is chosen here. 

In the following, we shall make use of eq. ( |C.8| ), and expand all the functions of z p in eqs. (|C.2|) 
and ( |U.3j) around 

2 



m B - m D 
m T /2 



0.0863, 



(C.12) 



neglecting C(e 3 ), 0(aj), 0(K i QCD /rrVQ) and 0(e x k 2 QCD /m 2 Q ). As argued in ref. J21], the leading 
renormalons cancel when such an expansion is performed for physical observables. Indeed, the 
QCD perturbation series converges remarkably well in our final result for the ratio (|C.1J): 



C 



g(z )\l + e 



7T 



+ K 



[a s C F f 



+ 



^ (^ 2) (-o) -P^izo) + P^(z ) 2 -pV{zo)pV>) + K 



K,Xi 



+ 



24A 5 



m 



T 



(1 ~ ZqY 



kA 2 



0.534 



1 + 0.019e + (0.001 - 0.005/ M (^ ) ± 0.002)e 2 



[a s C F ) 2 a s 
8tt 



Ai 



Y 



2ti(z )' 



A, 



7.1 GeV 12 GeV 



0.575 x (1 ± 0.01p ert ± 0.02 Al ± 0.02 A ) 



where 



2^(1- v^) 



g'(z ) 



2.94. 



7.6 GeV 2 _ 
(C.13) 

(C.14) 



In the middle step, the error in the 0(e 2 ) term originates only from pffl(zo) = —12.6 ± 0.4 (see 
eq. ( p.5| )). In the last step, f u (zo) is set to —0.9, and an overall ±1% perturbative uncertainty 
is assumed, resulting from the fact that our estimate of f u {zo) is very rough, from the error in 
PuK z o) an d from missing perturbative higher orders. As far as the non-perturbative parameters 
are concerned, we use Ai = (-0.27 ± 0.10 ± 0.04) GeV 2 0,0 A 2 = 0.12 GeV 2 and eq. Q. 

15 It is consistent with Ai = (-0.196±0.065±0.072) GeV 2 recently extracted by CLEO from the B -> X s -f 
spectrum. 
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A delicate point in our calculation of C is the fact that the unknown 0{h? QCD /m 2 D ) corrections 
have been neglected in the numerator of eq. (|C.8| ). Their potential effect can be studied by 
replacing Ai — > A x ± 0(Aq GD ) /tUd in eq. ( p.!3|) . However, if such corrections were sizeable, 
they would affect the determination of Ai from the semileptonic spectrum in ref. |21[| . In the 
present paper, we assume that those corrections are included in the error of Ai (and in its 
central value). A further study of this point is warranted. 

It is worth mentioning that our final central value C = 0.575 can be reproduced from the 
ratio of eqs. (|C.2|) and (|C.3|) expanded in a s and Aj, for z p = 0.285 2 at NLO, and z p = 0.298 2 at 



NNLO. Those two values are consistent with z p = (0.29 ± 0.02) 2 which has been used in many 
previous analyses of B — > X s j. On the other hand, when z p = (0.29 ± 0.02) 2 is used in such a 
ratio, one obtains at NNLO 

Cdirect = 0.595 x (1 ± 0.08 ± 0.05), (C.15) 

where the first error comes from z p , and the second one from poor convergence of the pertur- 
bation series. The uncertainties here are much larger than in the result (|C.13| ) obtained with 
the help of the T expansion. 



Appendix D 



In this appendix, we present the analytical formulae and calculate the numerical values of 
the functions a(z) and b(z) that occur in the expression for K c (|3.7|). Those two functions of 
z = [m c jm^) 2 originate from the two-loop b — > S7 matrix elements of the 4-quark operators 
Pi and P2 ( |2.1| ). The relevant diagrams are shown in fig. |6|. They were calculated in ref. ||1 1|| . 
Their colour structure implies that (Pi) = — g(i"2)- Additive constants in the functions a(z) 
and b(z) are chosen in such a way that both functions vanish at z — 0. Explicitly, 



a(z) 



b(z) 



f{ 

+ 



1^2 
3" 



3C(3) + (I 



fvr 2 L + \L 2 + ^ 



4 ~ 3" 



\, 2 L 



1 7"2 1 J_ r3 
4 -L + 12 L 



z 2 + 



2L - \L 2 



2 - |vr 2 + \L + \L 2 ) z+ (l-l^-L+ \L 2 ) z 2 + z 3 } } + O(z^), (D.l) 

{ (-3 + |tt 2 - L) z + (i + vr 2 - 2L - \L 2 ) z 2 + (-f 2 - |tt 2 - §L + 2L 2 ) z* 
-IttV/ 2 + m [-z + (1 - 2Z> 2 + (-f + f L) } + 0(z 4 L 4 ), (D.2) 



where L = \nz. 

The imaginary parts of a(z) and b(z) have very little influence on our prediction for BRy. 
Since the LO amplitude is real, they affect the r.h.s. of eq. (|3.5| ) only via 0(a 2 ) terms that 
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we neglect anyway^] and via the very small 0(V u b) correction in eq. ( |3.7| ). If we included the 
imaginary parts of K c and K t , our final result (|4.14j) would get enhanced by only around 0.5%. 




b s b s b s b 

Figure 6: Leading contributions to the b 



57 matrix elements of P 1 and P 2 . 



As we have mentioned in the introduction, the choice of renormalization scheme for m c and 
rrib in (Pip) is very important for BR 7 . In principle, this choice is a NNLO issue that can be 
resolved only after calculating three-loop corrections to the diagrams of fig. ||. Since calculating 
finite parts of such diagrams would be a very difficult task at present, we have to guess what 
the optimal choice of m c and rrib is, on the basis of our experience from other calculations. 

All the factors of m c in (Pip) originate from explicit mass factors in the charm-quark propa- 
gators. In the real parts of (Pip), those charm quarks are dominantly off-shell, with momentum 
scale [i set by rrib- Actually, we are not able to decide whether this scale is rrib, \ m b or \rrib- 
Therefore, we shall vary \i between m c ~ ^rrib and mi,, and use m^ s (/i) in a(z) and b(z). 

As far as the factors of nib are concerned, they originate either from the overall momentum 
release in b — > S7 or from the explicit appearance of rrib in the 6-quark propagators. In the 
first case, the appropriate choice of rrib is a low-virtuality mass. In the second case, there is 
no intuitive argument that could tell us whether m^ ole or rrib(mb) is preferred. However, we 
think that as long as the three-loop diagrams remain unknown, the best choice is to set all the 
factors of rrib equal to m\ s in a(z) and b(z). The distinction between m\ s and m^ ole as well as 
the infrared sensitivity of the pole mass can be ignored here, because the uncertainties due to 
m c {n) are very large. 

In determining the optimal renormalization scheme for rrib in the ratio rn 2 /rril, we have to 
take into account that the considered perturbative amplitude arises as a particular term in the 
operator product expansion for the inclusive decay of the B meson. Thus, the ratio ml/ml 
should be understood as m 2 /Q 2 , where Q 2 is the squared "mean" momentum of the 6-quark 

mr 



inside the 5-meson. Its value is given by one of the "kinetic" b-quark masses [21 



e.g. 



[m 



IS\2 



) 2 . Given the large error in m c (/i), it does not matter which definition of the bottom 



kinetic mass is chosen here. 

It remains to determine the numerical value of m^ IS (/i)/m^ s . From m c 



0.10) GeV |29| and the RGE for m c (/i), we find 



A1S 



is 



0.22 ±0.04, 



nT'c{ m c) — (1-25 ± 



(D.3) 



16 



Except for the 0(a 2 s ) contributions to the ratio r^o). 
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which implies 



a(z) = (0.97 ±0.25) + z(1.01 ± 0.15), (D.4) 
b(z) = (-0.04 ±0.01) + z(0.09±0.02), (D.5) 

for \i varying between m c and m^. The uncertainty in m c (m c ) is practically irrelevant here, when 
compared to the dominant error that originates from a variation of fi. We also note that our 
central value m c jm^ = 0.22, is very close to the RG-invariant MS ratio m c (//)/m&(^) = 0.215. 



Appendix E 



Here, we discuss the function B(Eq) in eq. ( |3.5|) that arises from b — > s^g and h — > s^qq with 
q — u,d, s: 



B(E ) 



(O)eff , 



7T 



2Eo 

m b 



(3qq(E ). 



(E.l) 



The functions <j>ij{5) originate from the gluon bremsstrahlung b — > s^g |13,|14], 47|. For z, j G 
{1, 2, 7, 8} their explicit form is as follows: 



> 22 



^(5) 
^78 (5) 



16z 
~2T 
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±4Li 2 (l -5) 5(2 + 5) ]n<J + 81n(l-<y) - -5 3 + 35 2 + 75 

3 3 



(E.2) 
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(E.4) 
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(E.6) 
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; 11 — — V22, 

36 
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<>22, 



'17 
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T 
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1 



■>27, 



'18 



'27, 



'2x 



'27- 



(E.7) 



They are identical to the fij(5) used in many previous analyses of B — > X s 7, except for the 



case i = j = 7. The difference between (f)jj(5) and f77(5) from eq. (37) of ref. |L2| is given 
by the first three terms in eq. ( |E.4| ): the Sudakov logarithms (In 2 5, In 5) and a constant term 
that makes $77(8) vanish at 5 — > 1. For simplicity, we refrain from resumming the Sudakov 
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logarithms here, because we are interested only in energy cut-offs Eq < 2.1 GeV (5 > 0.1), for 
which |077(5)| < 1, i.e. the logarithmic divergence of (firi(5) at 5 — > is not yet relevant. 
The function G(t) that appears in eqs. ( |E.2| ) and (|E.3|) reads 



G(t) 



-2arctan 2 Jt/(A-t), 



for t < 4 



-7T 



! /2 + 2 ln 2 [(^ + Vt^4)/2] - 2m ln[(V* + v / i 77 4)/2], for t > 4. 



(E.8) 



In our numerical analysis, the parameter z = m 2 c /m1 entering those equations is set equal to 
(0.22±0.04) 2 , as determined in appendix D. In eq. ( |E.6|) , we follow ref. |l2j and use m b /m s = 50. 

The functions 0y with 3 < (i or j) < 6 have only a 0.1% effect on BR[I? — > X/y]^ >i. 6 GeV- 
Therefore, we shall not give them explicitly here. They can be read from the results of ref. Jl4 . 

The coefficients C^ cS (fib) are equal to those in eqs. (20) and (22) of ref. Hl2] . They read 
(C< = Cf for % = 1,...,6) 
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- 0.0380 rf 6 


- 0.0185 rf 7 
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-lF t (r\ 4- 64017 
2 -t w -r 121Q12 



-0.9135 rf 5 + 0.0873 rf 6 - 0.0571 rf 7 + 0.0209 rf 8 , (E.9) 



where x = m t (m t ) 2 / . The powers aj can be found in table [I], in section |3|. 

It is important to mention that the sensitivity of B(Eq) to m b via the argument of (f)^ 
is rather weak when E = 1.6 GeV. When m b is varied between 4.5 and 4.9 GeV here, our 
final result for BR[_B — > Xg^E^isGeV is affected by only around 0.7%. In section |, we have 



m 



IS ; 



in eq. flElQ 



neglected this uncertainty, and used - 

The term denoted by j3 qq -(E Q ) in eq. ( |E.1| ) stands for the contributions from b — > sqq-y 
transitions with q = u,d, s. Perturbatively, such effects on T[b — > X s ^] are suppressed by either 



a 



3,.. .AM 



< 0.05 



or 



v: s v ub ci 0) M 



V t * s V tb C^ eS (fi b 



0.002 



(E.10) 



with respect to the leading terms. Further suppression occurs when we restrict ourselves to high- 
energy photons ||47]. In our numerical analysis we have set f3 q q(Eo) to zero without including 
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any additional uncertainty, which we expect to be acceptable at present for an energy cut-off 
E > 1.6 GeV. 
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